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The material KNi2Se2 has recently been shown to posses a number of striking physical properties, 
many of which are apparently related to the mixed valency of this system, in which there is on 
average one quasi-localized electron per every two Ni sites. Remarkably, the material exhibits a 
charge density wave (CDW) phase that disappears upon cooling, giving way to a low-temperature 
coherent phase characterized by an enhanced electron mass, reduced resistivity, and an enlarged 
unit cell free of structural distortion. Starting from an extended periodic Anderson model and 
using the slave-boson formulation, we develop a model for this system and study its properties 
within mean-field theory. We find a reentrant first-order transition from a CDW phase, in which 
the localized moments form singlet dimers, to a heavy Fermi liquid phase as temperature is lowered. 
The magnetic susceptibility is Pauli-like in both the high- and low-temperature regions, illustrating 
the lack of a single-ion Kondo regime such as that usually found in heavy-fermion materials. 



Heavy-fermion materials exhibit a host of fascinating 
collective quantum behaviors, which have made them a 
major focus of ongoing research for over three decades 
[212]. The "standard model" of heavy-fermion behavior, 
as depicted in the famed Doniach diagram [3J, features 
competition between the heavy Fermi liquid phase, stabi- 
lized by the hybridization between conduction electrons 
and localized moments, and a magnetic phase, typically 
stabilized by the RKKY interaction between the local- 
ized moments. An intriguing possibility that has received 
comparatively little attention is the existence of charge 
order, rather than the usual magnetic order, in proximity 
to the heavy-fermion state. Mixed- valency systems [4], 
which contain a variable number of localized electrons 
per atomic site, are a natural place to look for such com- 
peting effects. If the fractional filling takes a commen- 
surate value, then Coulomb repulsion between electrons 
on nearby sites may induce charge density wave (CDW) 
ordering. Mixed valency has been studied recently in 
/-electron materials exhibiting heavy-fermion behavior 
[5HZ], as well as in the context of the related "charge 
Kondo effect" [SH]. However, in both of those cases the 
emphasis has generally been on the single-site valency 
as hybridization or interaction between electrons is in- 
creased, rather than on the possibility of collective CDW 
formation and competition of this charge order with the 
heavy fermion phase. 

The material KNi2Se2 has recently been shown to ex- 
hibit several remarkable physical properties |10j . many 
of which appear to be related to its mixed-valent na- 
ture. At high temperatures the material has high resis- 
tivity; the magnetic susceptibility is constant, indicating 
Pauli paramagnetic response; and structural analysis re- 
veals that the material has at least three distinct sub- 
populations of Ni-Ni bond lengths. Upon cooling below 
Tcoh ~ 20K, the resistivity rapidly decreases, the struc- 
tural distortions disappear, and the material enters a co- 



herent heavy-fermion state with effective electron mass 
m* ~ 10mo, eventually giving way to superconductivity 
below T c « 1 K. This material is also unusual in that 
an applied magnetic field induces virtually no response 
in the measured specific heat and magnetoresistivity, in- 
dicating that the low-temperature coherent phase does 
not arise from competition with local magnetic order as 
in typical heavy-fermion materials. Rather, it was pro- 
posed that the coherent state competes with a charge- 
nuctuating state, facilitated by the mixed valency of the 
Ni ions in KNi2Se2 10J. The proximity of this system 
to charge (rather than magnetic) order, the disappear- 
ance of the CDW as temperature is lowered, and the 
lack of common signatures of Kondo-like behavior in the 
magnetic response together distinguish this material from 
conventional heavy-electron systems, making it a new 
and fascinating example of the diversity of effects that 
can be observed in these materials. 

In this study, we present a theory that captures the key 
ingredients that characterize this system. At high tem- 
peratures, the quasi-localized electrons in our model form 
a CDW and pair with one another into singlet dimers, 
which explains the observed structural distortion and in- 
scnsitivity to applied magnetic field. As temperature is 
lowered, the CDW dissolves in a first-order transition di- 
rectly into a spatially uniform, correlated heavy-fermion 
state, without the intermediate single-ion Kondo regime 
that is usually observed in heavy-fermion materials. The 
details of this model and the main results of the calcula- 
tions are presented below. 

KNi 2 Se2 has a quasi-two-dimensional, layered struc- 
ture, with the Ni and Se ions alternating in checkerboard 
fashion on a square lattice within each layer. Considera- 
tion of the stoichiometry reveals that the effective valency 
of Ni in this compound is "1.5+," so that at low ener- 
gies the effective degree of freedom is one quasi-localized 
d-electron with spin 1/2 per every two Ni sites, with a 
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small amplitude for these electrons to hop to neighboring 
Ni sites. Conduction electron bands are formed from the 
other Ni and Se orbitals and have a significantly greater 
bandwidth than the quasi-localized d-electrons. 

With this picture in mind, the following Hamiltonian 
describing the "extended periodic Anderson model" pro- 
vides a useful starting point: 

#ea = - t c + H - c -) - tfJ2( d l*dj„ + H.c.) 

(ij) (ij) 

-s f J2 n di + V X!(4x c «x + H.c.) ( X ) 
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+ U ndifridii + ^2 W tJ n di n d: j, 

where i,j denote Ni sites on a two-dimensional square 
lattice. The first term in this equation describes hopping 
of the conduction electrons. The second and third terms 
describe the hopping and on-site energy of quasi- localized 
electrons on neighboring Ni sites. The fourth term de- 
scribes hybridization between the two types of electrons. 
Finally, the last two terms describe Coulomb repulsion 
of d-electrons occupying the same site and nearby sites, 
where ndia = d icr di a . Summation over repeated spin in- 
dices is implicit in ([I]). The Hamiltonian is identical to 
the well-known periodic Anderson model, with the addi- 
tion of the W term describing intersite Coulomb repul- 
sion. This term is typically neglected in describing heavy 
fermion materials due to the fact that such systems usu- 
ally have exactly one local moment per site, so such a 
Coulomb term effectively adds an overall constant to the 
total energy. Such a term is crucial for a system such as 
KNi2Se2 near one quarter filling, however, since such sys- 
tems are susceptible to Coulomb repulsion-driven charge 
ordering. 

In the limit of large on-site repulsion U, it is convenient 
to enforce the constraint of no double occupancy through 
the introduction of slave boson operators [LTj - TTS] . In this 
formulation, we substitute di a = b\fi a , where fi a de- 
scribes a charge-neutral "spinon" that carries the spin of 
the electron, and the slave boson operator bi describes a 
spinless particle with positive charge. Within this formal- 
ism, rather than having to enforce the constraint < 1 
at every site, one instead has the constraint rifi+b\bi = 1, 
where n/, = f^ftf + f}±fn is the number operator for 
spinons. In terms of these new operators, the Hamilto- 
nian ([I]) becomes 



H = H C + H fc + H f + H w + H x + Hj (2) 
H c = -t c J2(cic^+H.c.) 
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FIG. 1. (a) T = diagram of CDW phases for various val- 
ues of second- and third-neighbor Coulomb repulsion, with 
J = V — 0. (b) Schematic illustration of the "plaquette" 
phase showing expected lattice distortion once the coupling 
between electron and lattice degrees of freedom is taken into 
account. Links connecting sites containing localized electrons 
are shorter and are represented by thicker lines, (c) Schematic 
illustration of a possible configuration in the case where long- 
range CDW order is absent, but the singlet nature of the 
state and distribution of bond lengths remain similar to those 
shown in (b). 



H w =^2W ij nf i nf j 
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The first four terms in |2]) are analogous to terms appear- 
ing in ([!]), but rewritten in the slave boson description. 
H\ replaces the on-site repulsion term in ([T]) by enforc- 
ing the constraint of no double occupancy for spinons 
via the Lagrange multiplier field A;. The last term in 
([2]) describes an antiferromagnetic Heisenberg interaction 
driven by superexchange between spins on neighboring 
sites. The Hamiltonian ^ is identical to the "Anderson- 
Heisenberg" model that has been studied recently [TH - 
[T5] . with the addition of the Coulomb term Hw- 

Figure [lja) shows the phase diagram with the pos- 
sible CDW phases for various values of first-, second-, 
and third-nearest neighbor repulsion (Wi, W2, and W3, 
respectively). We choose parameters such that the "pla- 
quette" CDW phase in the upper part of Figure [l|a) 
is realized, since this phase naturally allows for dimer 
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formation between neighboring spinons. Note that, due 
to the frustration between the competing tendencies for 
large W\ and W%, the realization of this phase does not 
require especially large values of third-neighbor repulsion 
W3. For example, the plaquette CDW phase is realized 
in the idealized case of unscreened ~ 1/r Coulomb re- 
pulsion, for which the ratios of the repulsion terms are 
W 2 /Wi = 1/V2 and W 3 /Wx = 1/2. When coupling be- 
tween electronic and lattice degrees of freedom is taken 
into account, such a picture can also explain the distinct 
peaks in the distribution of bond lengths observed at 
T > Tech via neutron pair distribution function analysis 
[10] , since links containing a singlet pair can be expected 
to be shorter than other links. It can be seen from Fig- 
ure[TJb) that each unit cell contains 2 short bonds, 2 long 
bonds, and 4 bonds of medium length. This is consistent 
with the three peaks in the distribution of bond lengths 
observed in experiment, with the central peak larger than 
the others. While there is no clear experimental evidence 
of long-range spatial order such as that described here, 
we expect that the key features of this model — spatially 
modulated electron density and singlet formation at high 
temperatures, giving way to a spatially uniform coherent 
state at low temperatures — will remain valid even in the 
absence of long-range order, as illustrated schematically 
in Figure |T](c). 

We proceed to study the Hamiltonian ^ within mean- 
field theory. Denoting as sublattice A (B) the sites shown 
as (un)occupied in the figure, the average occupation 
number is taken to be (n/j) = nj + Ci^' where n/ is 
the average density of spinons per site, A is the CDW 
order parameter, and Q = ±1 on sublattice A (B). At 
mean-field level, the nearest neighbor and second neigh- 
bor Coulomb terms W\$ merely shift the chemical po- 
tential for the spinons, so they will not be considered 
further here. The Lagrange multiplier field and the 
slave boson field bi are also treated as staggered mean 
fields: i\i — \q + QXi and &, = bg — Ci^i, with Ao,i and 
bo ,i real. The Heisenberg spins Si are expressed in terms 
of spinons: Sf = \f ia v'^{ifip- l where a a are the Pauli ma- 
trices. Defining the mean fields Xa,a<,b,ab = (fLfja), 
with xab defined on links between sites on different sub- 
lattices, xb on links between two sites on sublattice B, 
and xa(A') 011 links between two sites on sublattice A 
in the x(y)-direction. It is found that the free energy is 
always lowered in the CDW phase by having only one of 
XAi Xa' nonzero, so that the four spins on each plaquette 
form two dimers. In the uniform phase, Xij takes the 
same value on all links. 

In solving the mean-field equations, we require that 
the density of spinons is fixed to rif — 0.5 spinons per 
site in the limit V = 0, which is accomplished by ap- 
propriately setting the on-site energy of the spinons e / . 
Once the hybridization V is nonzero, £/ remains fixed 
to this value, and in general n/ 7^ 0.5 once V is finite. 
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FIG. 2. Upper plots show the T — phase diagrams for V 
vs. W3 with J = Stf (a), and for J vs. W3 with V = 3tf 

(b) . Lower plots show the reentrant first-order transition to 
the charge-ordered phase as a function of W3 with V — 3tf 

(c) , and as a function of V with Wz = 1.754/ (d), and with 
J = 8tf in both plots. 



The chemical potential is set to keep the total density 
of particles in the system fixed at n c + nf = 1.3 per 
site, which remains fixed even for V ^ 0. The results 
that follow are not particularly sensitive to the choice of 
n c , so long as n c > nj, so that there are enough con- 
duction electrons to screen all of the local moments in 
the coherent phase. The ratio of hopping amplitudes 
for spinons and c-electrons has been set to tf/t c = 0.2, 
which is substantially larger than that found in typical 
heavy-fermion systems. This is a reflection of the fact 
that the localized moments in KNi 2 Se2 are d-electrons, 
which are less tightly bound to their atomic cores than 
the /-electrons that constitute the local moments in most 
other heavy-fermion systems. This relatively small ratio 
of bandwidths is also the reason for the rather modest 
effective mass enhancement of m* ~ 10mo |10j . which 
is 10 ~ 100 times smaller than that typically found in 
/-electron heavy fermions. 

The phase diagrams shown in Figure|2]illustrate the ex- 
istence of a reentrant, first-order transition from a CDW 
to a spatially uniform phase upon cooling. This behav- 
ior is rather unusual, as in most systems the phase that 
breaks translational symmetry is the ground state that 
is realized as T — > 0. While the reentrance occurs along 
the entire critical line in Figure [2ja) , it emerges along 
the critical line in Figure [2jb) only for J > 3i/ , growing 
in extent as J is increased. These values of J are rather 
large to be generated by superexchange alone, for which 
one expects J ~ t 2 JU. It has been suggested that addi- 
tional contributions might arise in similar contexts from 
other superexchange processes in the CDW phase [17] . 
or from RKKY interactions at low temperatures [TH1 HH] • 
Reentrant behavior reminiscent of that shown here has 
been seen previously in a theory of simple checkerboard 
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FIG. 3. Densities of states for the CDW phase at T = 0.15t/ 
(a) and the coherent phase at T — O.Olt/ (b), with W3 = 
1.75i/, J = 8tf and V = 3t/. The Fermi level is at e = 0. 
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FIG. 4. Uniform magnetic susceptibility with W3 = 1.75i/ 
and J — 8tf. The light and dark dashed lines show x(T) in 
the uniform phase (V = 44/ ) and in the CDW phase (V = 
2tf), respectively. The solid line shows x(T) for V — Stf, 
which exhibits a phase transition at T = T co h = 0.12t/. All 
plots are normalized to x(0) for the solid curve. 



CDW ordering at 1/4 filling in layered molecular crys- 
tals [19] , although in that case only a second-order tran- 
sition was found. The model presented here also exhibits 
a second-order transition, but only at higher tempera- 
tures than those shown in Figure [2j The mean field \A 
is nonzero throughout the CDW phase shown in Figure 
[5J up to T ~ 0.6i/, indicating that singlet formation is 
taking place between the localized spinons. 

For the parameters given above, there is a jump in 
the average spinon occupation per site n/ at the first- 
order transition. At W3 — 1.75t/, the occupation jumps 
from n f = 0.45 for T < T coh = 0.12*/ to n f = 0.53 for 
T > T co h. The reduced valency nj < 0.5 in the uniform 
phase is expected, indicating that some of the spinons 
become delocalized and join the sea of conduction elec- 
trons in forming a heavy Fermi liquid. The increased 
valency in the CDW phase is less expected, but is consis- 
tent with the lack of long-range CDW order observed in 
experiment [TU], since a long-range ordered state would 
be impossible at incommensurate filling. The valency of 
localized electrons, which is given by n/ on account of the 
constraint (n/i) + bf = 1, can be measured using tech- 
niques such as resonant inelastic X-ray scattering [20] . It 
would be interesting to test whether the predicted jump 
in nt at T = T co h could be observed experimentally. 

The densities of states in the two phases are shown 
in Figure [3] In the CDW phase, the well-defined peaks 
above and below the Fermi level clearly show that the 
spinon excitations are gapped. In contrast, the spinons 
contribute to the hybridization peak at the Fermi level in 
the low-temperature phase, as is typical in heavy-fermion 
materials. The relative magnitude of the peak in this 
case, however, is substantially smaller than that in /- 
electron materials. Comparing the values of the densities 
of states at the Fermi level in the two different phases, one 
finds an enhancement of ~ 3.0 in the normal phase rela- 



tive to the CDW phase. This can be compared with mea- 
surements on KNi 2 Se2, where an enhancement of ~ 3.1 
was observed in the electronic specific heat coefficient 7 
for T < T coh [TO]. 

The uniform magnetic susceptibility can be calculated 
as the derivative of magnetization with respect to ap- 
plied magnetic field. Figure [4] shows the susceptibility 
as a function of temperature, assuming that the Lande 
g-factors for conduction electrons and spinons are equal. 
The approximately constant susceptibility in the region 
T < 0A2tf = T co h corresponds to the Pauli suscepti- 
bility of the heavy Fermi liquid, in which the conduc- 
tion electrons are hybridized with the localized spinons, 
leading to an enhanced density of states near the Fermi 
level. At T co h, there is a jump in the susceptibility to 
a much smaller constant value, indicating that only the 
conduction electrons contribute to the susceptibility at 
T > T co h, while the spinons form singlet pairs. This 
is in contrast to the Curie susceptibility x ~ 1 /T that is 
typically observed at high temperatures in heavy-fermion 
materials. 

Rather than exhibiting a sharp step, however, the ex- 
perimentally measured x(T) remains approximately con- 
stant through the CDW transition [TU] . One possible ex- 
planation for this discrepancy is the Van Vleck contribu- 
tion to x(T'), which has not been included in our model. 
The possibility of a large contribution of this type in 
heavy-fermion materials has been considered previously 
|21) . Others have since investigated the Van Vleck con- 
tribution to the susceptibility and have found that, when 
multiple localized bands are approximately degenerate, 
one generally has \v ~ XPauii [2?hf2~T|. Thus if an un- 
hybridized band that was not included in our model has 
large spectral weight near one of the peaks in Figure |3][a), 
an increased Xv could compensate for the decrease in 
XPauii at higher temperatures, with the sum of the two 
terms remaining roughly constant. Calculating the pre- 
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cise Van Vleck contribution to the susceptibility would 
require a more detailed knowledge of the band structure, 
however, and so we leave this as an open question to be 
addressed in future work. 



In conclusion, we have provided a theoretical frame- 
work for describing the key properties of the recently 
discovered mixed- valency material KNi2Se2, most impor- 
tantly the vanishing of the CDW phase upon cooling. 
The formation of singlet dimers by the local moments 
in the CDW phase explains the lack of common signa- 
tures of single-ion Kondo behavior, such as a Curie sus- 
ceptibility at high temperatures. This mechanism may 
also explain the lack of a resistivity peak in measure- 
ments on KNi 2 Se2 [ID] . Such a peak typically forms 
in heavy-fermion materials at temperatures just above 
the coherence temperature, where the Kondo screening 
clouds act as spin-flip scattering centers but are not yet 
coherent with one another. A direct transition from a 
singlet CDW phase to a coherent low-temperature phase 
precludes this possibility however, and is consistent with 
the monotonically decreasing resistivity observed in ex- 
periment as T is lowered. This material illustrates the 
potential of mixed-valency systems for exhibiting a rich 
array of collective quantum behaviors. 
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